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Abstract 

We report measurements of the differential q 2 spectrum and the forward-backward asymmetry 
for B — ► K^l + l~ , where I represents an electron or a muon, with a data sample of 253 fb -1 accu- 
mulated on the T(45) resonance with the Belle detector at KEKB. We also present measurements 
of the branching fractions and their ratios. 

PACS numbers: 13.20.He, 11.30.Hv 



3 



Flavor-changing neutral current (FCNC) processes are forbidden at tree level in the Stan- 
dard Model (SM); rather, they proceed at a low rate via loop or box diagrams. If additional 
like diagrams with non-SM particles contribute to such a decay, their amplitudes will inter- 
fere with the SM amplitudes and thereby modify the decay rate as well as other properties. 
This feature makes FCNC processes an ideal place to search for new physics. 

Measurements of the radiative penguin decay B — > X s j 0, 0, 0], which are consistent 
with the SM prediction, strongly constrain the magnitude — but not the sign — of the effective 
Wilson coefficient Cj. This is an important limitation, since non-SM contributions can 
change the sign of C-j without changing the B — > X s 7 branching fraction 0]. 

The b — > s£ + £~ process is promising from this point of view, since not only the photonic 
penguin diagram but also the Z-penguin and box diagrams contribute to this decay mode. 
As a result, the Wilson coefficients C 7 , C 9 and C w can be completely determined. The first 
observations of B — > K£ + £~B — > K*£ + £~ and inclusive B —>■ X s £ + £~ decays were reported 
by the Belle Collaboration jH S Q ■ The measured branching fractions of these decay modes 
were used to exclude a large area of the allowed region in the C 9 -C w plane (H 0| . However, 
the determination of the sign of CV (as well as of C 9 and Ci ) requires precise measurements 
of the distribution in squared dilepton momentum q 2 and the forward-backward asymmetry 
in these decay modes [10j|. 

In this paper, we present preliminary results of improved measurements of B — > K£ + £~ 
and B — > K*£ + £~ using data produced in e + e~ annihilation at the KEKB asymmetric 
collider [n] and collected with the Belle detector. The data sample corresponds to 253 fb _1 
taken at the T(4S) resonance and contains approximately 275 million BB pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere The data were collected with two 
different inner detector configurations. For the first sample of 152 million BB pairs, a 2.0 
cm radius beampipe and a 3-layer silicon vertex detector were used; for the latter 123 million 
BB pairs, a 1.5 cm radius beampipe, a 4-layer silicon detector, and a small-cell inner drift 
chamber were used |l3j |. 

In this analysis, primary charged tracks — except for the daughters from Kg — > 7i + tt~ 
decays — are required to have impact parameters relative to the interaction point of less 
than 5.0 cm along the z axis (aligned opposite the positron beam) and less than 0.5 cm 
in the r<p plane that is transverse to this axis. This requirement reduces the combinatorial 
background from photon conversion, beam-gas and beam-wall events. 

Charged tracks are identified as either kaons or pions by a likelihood ratio based on the 
CDC specific ionization, time-of-flight information and the light yield in the ACC. This 
classification is superseded for a track that is identified as an electron or for a pion-like track 
that is identified as a muon. Electrons are identified from the ratio of shower energy of the 
matching ECL cluster to the momentum measured by the CDC, the transverse shower shape 
of this cluster, the specific ionization in the CDC and the ACC response. We require that 
the electron momentum be greater than 0.4 GeV/c to reach the ECL. Muons are identified 
by their penetration depth and transverse scattering in the KLM. The muon momentum 
is required to exceed 0.7 GeV/c. The muon identification criteria are more stringent for 
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momenta below 1.0 GeV/c to suppress misidentifled hadrons. 

Photons are selected from isolated neutral clusters in the ECL with energy greater than 
50 MeV and a shape that is consistent with an electromagnetic shower. Neutral pion can- 
didates are reconstructed from pairs of photons, and are required to have an invariant mass 
within 10 MeV/c 2 of the nominal 7r° mass and a laboratory momentum greater than 0.1 
GeV / c. The pion momentum is recalculated by constraining the invariant mass to the nom- 
inal 7T° mass. Kg candidates are reconstructed from oppositely charged pions that have 
invariant masses within 15 MeV/c 2 of the nominal Kg mass. We impose additional cri- 
teria based on the radial impact parameters of the pions (Sr), the distance between the 
closest approaches of the pions along the beam direction (Sz), the distance of the vertex 
from the interaction point (/), and the azimuthal angle difference between the vertex di- 
rection and the K s momentum direction (5(f)). These variables are combined as follows: 
for p K o < 0.5 GeV/c, 5z < 8 mm, 5r > 0.5 mm, and 5<f> < 0.3 rad are required; for 
0.5 GeV/c < p K o < 1.5 GeV/c, 5z < 18 mm, 5r > 0.3 mm, 5(f) < 0.1 rad, and / > 0.8 mm 
are required; and for p K o > 1.5 GeV/c, 5z < 24 mm, 5r > 0.2 mm, 5(f> < 0.03 rad, and 
/ > 2.2 mm are required. 

K* candidates are formed by combining a kaon and a pion: K + n~ , K s ir + or K + n° 
The K* invariant mass is required to lie within 75 MeV/c 2 of the nominal K* mass. For 
modes involving neutral pions, combinatorial backgrounds are reduced by the additional 
requirement cos#hei < 0.8, where ^hei is defined as the angle between the opposite of B 
momentum and the kaon momentum direction in the K* rest frame. 

B candidates are reconstructed from a K^ candidate and an oppositely-charged lepton 
pair. We use two variables defined in the CM frame to select B candidates: the beam-energy 

constrained mass Mb c = \J E£ CSLIQ 2 — p* B 2 and the energy difference AE = E* B — E^ CSLIQ , where 
p* B and E B are the measured momentum and energy, respectively, of the B candidate, and 
-^beam * s the beam energy (T^J. When multiple candidates are found in an event, we select 
the candidate with the smallest value of |A£J|. 

Backgrounds from B — > J /if) (if)') K^*' are rejected using the dilepton invariant mass. The 
veto windows are defined as 

lj/i> 

-0.20 GeV/c 2 < M e+e - - My < 0.07 GeV/c 2 for K* and K modes 
-0.15 GeV/c 2 < M M+M - - Mj^ < 0.08 GeV/c 2 for K* modes 
-0.10 GeV/c 2 < M M+M - - Mj^ < 0.08 GeV/c 2 for K modes 
-0.10 GeV/c 2 < M M+M - - My < 0.08 GeV/c 2 for K* and K modes 
If a photon with energy less than 500 MeV is found in the 50 mrad cone along the electron mo- 
mentum direction, we reapply the above vetoes with the invariant mass calculated including 
this photon to reject the background with a bremsstrahlung photon, J/if)(ip') —>■ e + e~ r ]( r y). 
For K*e + e~ modes, B — > J/if) K can be a background if a bremsstrahlung photon is missed 
and a pion from the other B meson in the event is included. We suppress this background by 
the following prescription: the included pion is discarded and an unobserved bremsstrahlung 
photon is added with a direction parallel to the electron or positron and an energy that gives 
AE = for the B candidate. If the dilepton mass and the beam-energy constrained mass 
are consistent with a B — ► J /if) K event, the candidate is vetoed. 

We suppress background from photon conversions and ir° Dalitz decays by requiring the 
dielectron mass to satisfy M e + e - > 0.14 GeV/c 2 . This cut eliminates possible peaking 
background from B — > K*^ and B — > K^tt . 



-0.25 GeV/c 2 < M e+e - - M J/V , < 0.07 GeV/c 2 for K* modes 
-0.20 GeV/c 2 < M e+e - - M J/V , < 0.07 GeV/c 2 for K modes 
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Background from continuum qq events is suppressed using event topology. Continuum 
events have a jet-like shape while BB events have a spherical shape in the center-of-mass 
frame. A Fisher discriminant T |l6j is calculated from the energy flow in 9 cones along the 
B candidate sphericity axis and the normalized second Fox- Wolfram moment R2 (l7| - We 
combine this with the cosine of the polar angle cos 9* B of the B meson flight direction and the 
cosine of the polar angle cos 9* ph of the B meson sphericity axis to define likelihoods £ sig and 
£ CO nt f° r sig na l an d continuum background, respectively, and then cut on the likelihood ratio 
T^-cont = £ S ig/(£ S ig + ^cont)- For the muon mode, | cos#* ph | is not used since its distribution 
is nearly the same for signal and continuum within the detector acceptance. 

The dominant background from BB events is due to semileptonic B decays. The missing 
energy of the event, -E m i SS = 2E£ eaja — E* is where E* is is a total visible energy in the event, 
is used to suppress this background since the undetected neutrinos carry away a substantial 
amount of energy. The B meson flight angle, cos 9%, is also used to suppress combinatorial 
background in BB events. We combine -E m i SS and cos 6^ into signal and i?5-background 
likelihoods and cut on the likelihood ratio 1Z BB , defined similarly- to 7£ CO nt- 

The B — > K(*H + £~ decays are generated according to Ref [8| and then processed by a 
GEANT-based Monte Carlo (MC) simulation to estimate the efficiencies. The signal box 
is defined as |M bc - M B \ < 0.007 GeV/c 2 for both lepton modes and -0.055 GeV< AE < 
0.035 GeV (|A£7| < 0.035 GeV) for the electron (muon) mode. We make distinct selections 
on 7Z con t and 7Z BB for each decay mode. The detection efficiencies are estimated from the 
MC samples and are summarized in Table H] 

To determine the signal yield, we perform a binned maximum-likelihood fit to each 
distribution. The expected number of signal events is calculated as a function of M^ c using 
a Gaussian signal distribution plus background functions. The mean and the width of the 
signal Gaussian are determined using observed J/ipK^*' events. A MC study shows that 
the width has no dependence on the dilepton invariant mass. We consider backgrounds 
from three sources: semileptonic decays, decays containing J/ip and ip' mesons, and double 
misidentification of hadrons as leptons. The background from semileptonic decays is parame- 
terized by the ARGUS function |20j. The shape is determined from large BB and continuum 
MC samples, each containing at least one oppositely charged lepton pair. The shape param- 
eter obtained from MC is consistent with that taken from a data sample of B — > K^e ± fi T 
candidates. The residual background from J/ip and ip' mesons that cannot be removed by 
the tp^ veto windows is estimated from a large MC sample of J/ip and ip' inclusive events 
and parameterized by an ARGUS function and a Gaussian. The background contribution 
due to misidentification of hadrons as leptons is parameterized by another ARGUS function 
and a Gaussian. The ARGUS function represents the combinatorial background while the 
Gaussian represents the component that forms a peak in the M^ c distribution. The shape 
and normalization of this background are fixed using the B — > K^*'h + h" data sample (where 
h refers to a pion or kaon). All K^*'h + h~ combinations are weighted by the momentum- 
and polar angle- dependent probability of misidentifying K^*'h + h~ as K^*>£ + £~ . This study 
yields 1.06 ± 0.09 Kh + h~ events and 0.60 ± 0.05 K*h + h~ events in the peak region. Other 
backgrounds with misidentified leptons are negligible. The normalizations of the signal and 
the background from real leptons are floated in the fit. The fit results are shown in Fig. ^ 
and Table |U 



We observe 78.5l}g;J B K*£+£~ signal with a significance of 11.0 and &2.2±\l$ B 



K£ + i signal with a significance of 10.1. The significance is defined as J — 2 ln(£ /^max), 
where £ max is the maximum likelihood in the Mb c fit and Co is the likelihood of the best 
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M bc (GeV/c 2 ) 



FIG. 1: M bc distributions for (a) B -» Ke + e~, (b) 5 -> iTe+e", (c) B -» KfJL + fi~, (d) B -» 
K*[i + n~, (e) £? — > Kl + £~ and (f) S — > samples. The solid and dashed curves are the fit 

results of the total and background contributions. 
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fit when the signal yield is constrained to be zero. Since the signal and background shapes 
are fixed in the fitting, the error due to uncertainty in shape parameters is included in the 
systematic error. We determine this by varying each shape parameter by ±ler and recalcu- 
lating the signal yield. We quote the smallest resulting significance in our determinations 
above. 

TABLE I: Summary of the fit results and branching fractions. Signal yield estimated from the M^ c 
fit, detection efficiency for each mode, obtained branching fraction, 90% confidence level upper 
limit of the branching fraction and the significance of the signal. The first error in the signal 
yield and branching fraction is statistical, the second is systematic and the third is due to model 
dependence. The first error in the efficiency is due to MC statistics and systematic effects and the 
second one is due to model dependence. 



Mode Signal yield Efficiency [%] B [xl(T 7 ] U.L.[xl(r 7 ] Signif. 

K°e + e~ — 1.0±J;|±3:| 5.10 ± 0.29 ± 0.20 -0.70lJ;||±^| ± 0.03 3.0 0.0 

K + e + e- 28.6±f;J±g:f 16.3 ±0.7 ±0.1 6.40lJ;|glg;g ± 0.05 - 6.4 

Ke+e" 26.6±|;f±g;l 10.7 ±0.5 ±0.1 4.541^1^ ± 0.06 - 5.5 

K*°e + e~ 22.0±f;|±g;? 4.32 ± 0.22 ± 0.36 18.5±|;| ± 1.1 ± 1.5 - 4.4 

K* + e+e- 6.2±f;£±8;jj 1.42 ± 0.08 ± 0.06 16.0^ 7 4 ±L8 ± °- 7 37 L7 

K*e + e~ 28.9^j ±0.9 2.87 ±0.15 ±0.21 18.4±|;| ± 1.1 ± 1.3 - 4.8 

K°fi + fi- H.6±|;2±o:3 6.76 ±0.41 ±0.04 6^26±f3I±^|frbToTo4 - 5A~ 

K+n+Lt- 39.9±|;|±g;f 23.2 ±1.1 ±0.5 6.28i^±g;|? ± 0.13 - 8.3 

Kn + ir 51.5±f;|±g;f 15.0 ±0.7 ±0.2 6.26±5;g|±g;|^ db 0.10 - 9.7 

K*°n + n~ 40.7l^glo;6 8.01 ±0.43 ±0.14 18.5±|;f ± 1.0 ± 0.3 - 8.4 

K* + n+^- 11.4±f;|±g;§ 2.55 ±0.14 ±0.08 16.3±|;|±?;| ± 0.5 - 3.5 

IC/j+fT 52.5i|;Jlg;| 5.28 ±0.29 ±0.03 18.1±|;^ ± 1.1 ± 0.1 - 9.1 

K°£ + £~ 10.7±|;|±g;§ 5.95 ± 0.36 ± 0.11 3.28±};f|±g;|^ ± 0.06 - 

K+£+£- 68.6±|;|±?;g 19.8 ± 0.9 ± 0.2 6.32±g;|§ ± 0.29 ± 0.06 - 9.4 

K£+£- 78.5tjo:oii:! 13.0 ±0.6 ±0.4 5.50^;^ ± 0.27 ± 0.02 - 11.0 

K*°£+£- 63.8±|;| ±0.9 6.88 ±0.37 ±0.30 16.9±1;| ± °- 9 ± °- 7 " 9 - 4 

K*+£+£- 16.3±2;|±g;| 2.22 ± 0.13 ± 0.02 13.4±^;g±?;g ±0.1 3.7 

K*£+£~ 82.2±iJ:|±i;2 4.55 ± 0.24 ± 0.12 16.5±|;| ± 0.9 ± 0.4 - 10.1 



We consider experimental systematic effects from the fit, the efficiency determination 
and BB event counting. Uncertainty in the background function is the dominant source 
of the systematic error. To evaluate the change caused by the uncertainty in the signal 
function parameters, the mean and the width of the Gaussian are changed by ±lcr from 
the values determined from J/ipK^ events. The uncertainty in the background shape is 
obtained by varying the ARGUS shape parameter by ±lcr from the value determined with 
a large MC sample. The uncertainty in the peaking background contribution is evaluated 
by changing the area of the associated Gaussian by ±lcr. The systematic errors associated 
with the fit function are shown in the second column of Table |l] Systematic uncertainties in 
the tracking, charged kaon identification, charged pion identification, electron identification, 
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muon identification, Kg detection and 7r° detection efficiencies are estimated to be 1.0%, 
1.0%, 0.8%, 0.5%, 1.2%, 4.5% and 2.7% per particle, respectively. The uncertainty in 
the background suppression is estimated to be 2.3% using J/ipK^ control samples. The 
systematic error due to MC statistics is less than 0.7%. The uncertainty in BB event 
counting is 1.1%. The systematic errors associated with efficiency and BB event counting 
are summarized in Table |HI Total experimental systematic errors are calculated by adding 
all systematic errors in quadrature. 

The systematic uncertainty due to theoretical modeling is also evaluated. We apply the 
same selection criteria on the signal samples generated according to the three form factor 
models [H, EH, 0] and obtain the efficiencies. The maximum difference in these efficiencies 
is assigned as uncertainty in model dependence and is listed as the final value in column 
four of Table |TJ 



Source 


K° 


K+ 


K+tt- 




K+TT° 


Tracking 


2.0 


3.0 


4.0 


3.0 


3.0 


Kaon identification 




1.0 


1.0 




1.0 


Pion identification 






0.8 


0.8 




Lepton identification (e//x) 


1.0/2.4 


1.0/2.4 


1.0/2.4 


1.0/2.4 


1.0/2.4 


Kg detection 


4.5 






4.5 




7T° detection 










2.7 


BG suppression 


2.3 


2.3 


2.3 


2.3 


2.3 


BB event counting 


1.1 


1.1 


1.1 


1.1 


1.1 


total (e/fi) 


5.8/6.1 


4.2/4.7 


5.1/5.5 


6.7/6.8 


6.3/6.1 



TABLE II: Summary of systematic errors in efficiencies and BB event counting. 

In the calculation of the branching fraction, we assume equal production rates of charged 
and neutral B meson pairs from T(4S") as well as isospin invariance. When combining the 
K*e + e~ and K* /x + /i~ modes, we assume the ratio of branching fraction of K*n + ji~ to that 
of K*e + e~ to be 0.75 0]. The combined branching fraction corresponds to the muon mode. 
The branching fractions are found to be 

B(B -> K£+£-) = (5.50±g ;g ± 0.27 ± 0.02) x 10~ 7 
B{B -> K*£ + £-) = (16.5l|| ± 0.9 ± 0.4) x lO" 7 , 

where the first error is statistical, the second is systematic, and the third is due to model 



dependence. These values are consistent with the SM predictions |8j, [18J, [19l |21(. The 
branching fraction of B — > K*£ + £~ is slightly larger than our previous result [6j. Since we 
now use NNLO effective Wilson coefficients |8j, which gives larger CV and smaller Cg values 
than the older NLO calculation, the efficiency of B — * K*£ + £~ is smaller by about 12%. As 
a result, the branching fraction is now higher. (With the NLO effective Wilson coefficients, 
the branching fraction becomes 14.7 x 10~ 7 , which is consistent with our previous result.) 
The branching fractions of the other decay modes are listed in Table |l] 

For the modes with a significance of less than 3.0, we also set 90% confidence level (C.L.) 
upper limits for the branching fractions. The 90% C.L. upper limit yield N is defined as 



/{f C(n)dn = 0.9 / °° C{n)dn. The function C(n) is the likelihood with signal yield n, where 
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each signal and background shape parameter is modified up or down by its systematic error 
in the direction that increases the signal yield. The upper limits for the branching fraction 
are then calculated by reducing the efficiency by its systematic error. 

We calculate the ratios of branching fractions to muon and electron modes. The ratio 
of branching fraction of B — > Kfi + fi~ to B — > Ke + e~ (TZku)^ sensitive to neutral Higgs 
emission from the internal loop in the two Higgs doublet model with large tan (3 2^1 . If the 



Higgs contribution is sizable, this ratio is greater than unity. The corresponding ratio for 
B —* K*£ + £~ (7Zk*m) is sensitive to the size of photon pole and is predicted to be about 
0.75 in the SM. The ratios are measured to be 

T> — 1 qq+0.39+0.06 
K-KU — i -'JO-0.41-0.07 

n KHi = o.98±8:|S ± o.os, 

which are consistent with the SM predictions. 

The distribution of the squared dilepton momenum q 2 and forward-backward asymmetries 
in B — > K£ + £~ and K*£ + £~ are also measured. The forward-backward asymmetry (A^b) is 
defined as the partial rate asymmetry between the positive and negative regions of cos9 Bi +, 
the cosine of the angle between the B° or B + meson and positive charged lepton momentum 
directions in the dilepton rest frame, 



.4 



FB 



r(cos# B£+ > o) - r(cos0 B£ + < o) 
r(cos0 B £+ > o) + r(cos^+ < 0)' 



The signal yield is extracted from a fit to the M bc distributions in each q 2 bin and cos 6^+ 
region. This procedure takes into account the forward-backward asymmetry of the back- 
ground. The q 2 resolution is about 0.6% and this is small enough relative to the bin size. 
For the q 2 distribution measurement, the efficiency in each q 2 bin is obtained from the MC 
samples described earlier. The branching fraction at low q 2 differs between K*e + e~ and 
K*fi + fi~; we use the average in calculating the K*£ + £~ differential branching fraction. Fig- 
ure^ and Table UTTI show the q 2 distributions; they are in agreement with the SM predictions. 

Figure H3 shows the raw forward-backward asymmetry distributions for B — > K£ + £~ and 
K*£ + £~ . The asymmetry in B — > K£ + £~ is expected to vanish in the SM; this expectation 
is essentially unchanged by the presence of new physics 23]. Therefore, B — > K£ + £~ is a 
good control sample for the more interesting forward-backward asymmetry in B — > K*£ + £~ . 
The measured asymmetry for B — > K£ + £~ is indeed consistent with zero. For B — > K*£ + £~ , 
the asymmetry is in agreement with both the SM and the wrong-sign C-j expectations, since 
the statistical power is not yet sufficient to distinguish between these two cases. 

In summary, we report first measurement of forward-backward asymmetry as a function 
of q 2 in B — ► K*£ + £~ . Within the limited statistical precision, the measured asymmetry in 
B — > K*£ + £~ is consistent with both the SM and the wrong sign C7 case. We also report 
improved measurements of the branching fractions of B — > K£ + £~ and B — K*£ + £~ . The 
measured values, their ratios and the q 2 distributions for B — > K£ + £~ and B — > K*£ + £~ are 
in good agreement with the SM predictions [1, 0, 1^, 21 . 
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20 

q 2 (GeV 2 /c 2 ) 



20 



FIG. 2: The q 2 distributions of (a) K£ + £ and (b) K*£ + £ . Points with error bars show the data, 
while the hatched bands show the range of the SM predictions 0, H^ . Il9| . 



K£+£~ 








q 2 /GeV 2 


yield 


efficiency [%] 


B[xl(T 7 ] 



[0,4] 

[4,8] 

[8,12] 

[12,16] 

[16,20] 

[20, 24] 



+4.4 +0.3 
-3.7 -0.4 
-6.1 
-5.4 

+4.5 +0.4 
-3.8 -0.5 

H ±0.3 

2a4+5 .4+0.4 



12.4 
26.7: 
10.1 
6.5 1 



,+0.36 



-0.30 



±0.06 ±0.01 



-4.8 -0.5 

1.3 ±^±0.3 



11.3 ±0.5 ±0.1 1.00 
± 0.4 17 

12.0 ±0.6 ±0.1 

11.1 ±0.5 ±0.3 0.53 1^ ± 0.04 ± 0.01 
17.0 ±0.8 ±0.3 
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.0 ± 0.8 ± 0.3 1.43 toil ± 0.07 ± 0.01 

n 7fi +0.34 +0.05 +n m 
u - ' D -0.29 -0.06 ^ U - Ui 
, +0.30 
-0.24 

1-09 l° ; 2 2 9 6 i^± 0.02 
2 ± 0.6 ± 0.8 0.09 Iftio ± 0.02 ± 0.006 



K*£ + £- 








q 2 GeV 2 


yield 


efficiency [%] 


B[xl0~ 7 ] 



[0,4] 
[4,8] 
[8, 12] 
[12, 16] 
[16,20] 



11.3±f;| ±0.6 2. 

t-6.1 +0.4 a 
-5.4-0.6 



21.6; 



16 ± 0.11 ± 0.56 4.8 tH ± 0.4 ± 1.2 
09 ± 0.22 ± 0.11 4.8 t\i ± 0.3 ± 0.3 
83 ±0.15 ±0.18 5.0 jl^-J ± 0.2 ± 0.1 
71 ± 0.14 ± 0.07 4.3 ± 0.2 ± 0.1 
I G.O nx ± 0- 5 3 - 99 ± 0- 21 ± 0- 26 3 - 6 -\l ± 0- 2 ± 0- 2 



15.7i||±0.4 2. 
mill; 7 , ±0.4 2 ' 

,+5.1 



TABLE III: Yield, efficiency and branching fraction in each q 2 region. The first error is statistical, 
the second is systematic and the third is model dependence. 
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FIG. 3: The raw forward-backward asymmetry in (a) K£ + £~ and (b) K*£ + £~. Points with 
error bars show the data. Curves show theoretical predictions including the effect of experimental 
efficiency and resolution for the SM case (blue) and wrong-sign CV case (red). The pairs of dashed 
lines are the boundaries of the J ftp and ip' veto regions. 
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